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ABSTRACT 


A MW class free electron laser capable of delivering energy at the speed of light 
can improve ASCM defensive capabihty for Navy ships. Many design challenges must 
be overcome to make such a weapon possible. One such challenge is to maintain the 
power density on laser cavity mirrors at acceptable levels. The use of a short Rayleigh 
length to increase beam spot size at the mirror is studied as a possible solution to this 
problem. In this thesis, it is shown that by using a short Rayleigh length PEL, power 
densities at the mir rors are significantly reduced without causing a noticeable reduction in 
performance. 

For a short Rayleigh length FEL, the resonator cavity is sensitive to misahgnment 
and vibration. The effect of mirror tilt due to vibrations is explored and the results show 
that as mirror tilt increases, FEE efficiency does decreases. However, a mirror tilt several 
orders of magnitude greater than currently achievable active ahgnment tolerances is 
required before the FEE efficiency is noticeably affected. In this thesis, it is shown that 
mirror tilt within achievable tolerance limits will not adversely affect the performance of 
a FEE. 
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1. INTRODUCTION 


The U.S. Navy’s interest in a high energy laser was expressed by the Commander 
in Chief of the U.S. Atlantic Fleet on the 24^^ of April, 2001, in a letter to the Chief of 
Naval Operations: 

...that the Speed of Light weapons can be very effective against these 
small high speed threats. Such a laser weapon would offer our Naval 
forces an extremely versatile weapon to counter numerous soft and hard 
targets. A High Energy Laser weapon can be designed to dehver energy 
that can track, warn, damage, mission kill, and if need be, destroy a threat. 

I beheve it is exactly this type of weapon system that our forces need in 
the httoral environment where, even though the threat may not be as 
sophisticated as a highly maneuverable cmise missile.... 

The free electron laser (EEL) is a possible answer for the Navy’s desire for a high 
energy laser (HEL) weapon system. Chapter n of this thesis presents background 
information on why directed energy weapons, and specifically the EEL, are a possible 
choice as an anti-ship cmise missile defense system. 

Chapter in gives a detailed description of the Thomas Jefferson Nuclear 
Accelerator Eacihties (TJNAE) proposed 100 kW EEL. The description includes a brief 
overview of EEL components and the parameters used for this particular design. 

Chapter IV presents a basic overview of EEL theory and physics required to 
understand the operations of a EEL. The description of how energy is transferred from 
the electron beam to the laser beam is presented along with the supporting theory and 
equations. 

Chapter V presents simulation results of short Rayleigh length EEL’s. A short 
Rayleigh length EEL is explored as a possible method to decrease the power density on 
laser cavity optical mirrors, thus preventing mirror damage. Due to the large number of 
simulations required for this research, the work was divided among several people with 
my contribution consisting of simulations and analysis for Rayleigh lengths of Zq = 0 IL 
and 0.2L where L is the undulator length. This material was originally presented at the 
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23 International FEL Conference in Darmstadt Germany in August 2001 and has been 
published in Nuclear Instmments and Methods in Physics Research [1]. 

Chapter VI presents simulation results describing stabihty analysis of the FEL 
when operating with a slight mirror tilt. This is the first such study of the effect of mirror 
tilt on the performance of the FEL. For an FEL using a short Rayleigh length, the optical 
mode may be sensitive to mirror tilt angle due to vibrations. Multiple simulations were 
conducted to determine the actual affect on FEL performance. Again, due to the large 
number of simulations, the work was divided among several people, with my contribution 
consisting of simulation and analysis of data for mir ror tilt angles from 0 to 400 irad. 
This material was originally presented at the 24'*’ International FEL Conference in 
Chicago, Ill inois in September 2002, and will be pubhshed in 2003 

Chapter VII presents an analysis to determine the required power for a shipboard 
point defense system FEL. Target kill mechanisms, and atmospheric propagation losses 
are considered to show that a MW class FEL has the potential to be a successful ASCM 
point defense system. 


2 



n. BACKGROUND 


A. WHY DIRECTED ENERGY 

Since the end of the Cold War, the U.S. Navy’s focus has shifted from preparing 
for warfare on the open ocean to developing capabihties for conducting combat 
operations in the littoral waters of the world. The prohferation of increasingly 
sophisticated anti-ship cmise missiles (ASCM’s) threatens the abihty of Navy ships to 
operate and survive close to hostile shores. These missiles are low cost weapons that can 
i nfli ct serious damage to a ship, and must be accounted for in the ship’s defensive 
capabihties. 

Current anti-ship missile defense systems include the Phalanx close-in weapons 
system (CfWS), RIM-116A roUing airframe missile (RAM), and standard missile. These 
systems provide some measure of defense for the ship against ASCM’s, but they have 
limitations that make ships vulnerable to anti-ship missile attacks. These systems are 
based on gunnery and rocketry, and have a time of flight required to get to the incoming 
threat and destroy it. There are a lim ited number of weapons available in the magazine 
for use, and once they are depleted the ship is defenseless. The lim it of defensive 
capabihties with these types of defensive weapons is rapidly being approached while the 
technology for ASCM’s has continued to advance. Short range detection of high-speed 
threat missiles may result in the missiles being destroyed in close proximity to the ship. 
Destroying the missile minimizes damage to the ship, but missile fragments may strike 
the ship and cause significant damage to unprotected equipment and personnel. Figure 1 
shows the probabihty of missile fragments hitting the ship as a function of range. For 
this simulation, a typical missile travehng at 400 m/s at an altitude of 8 m with a spread 
of missile fragment sizes based on known missile break up from weapons testing is used. 
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Probability of Missile Fragment Hit vs Range 



Figure 1. Probability of Missile Fragment Hitting Ship 

ASCM’s have become faster and more maneuverable, greatly increasing the risk 
to US warships. The Russian SS-N-26 Oniks is a Mach 3.5 sea skimming ASCM and is 
an example of velocities representative of these new threat missiles [2]. As the speed of 
the threat missile increases, the reaction time a ship has to identify an incoming threat, 
track, develop a solution, and launch a defensive weapon gets shorter. These faster 
ASCM’s and shorter detection ranges are stressing current defense systems to their lim it 
and the only way to keep up with the threat is to pursue new technologies that offer the 
advantages of faster speed. One such technology currently being considered is directed 
energy weapons. Directed energy in the form of a laser offers “speed of tight” capability 
and the ability to put large amounts of power onto an incoming missile in a short amount 

of time, destroying the target at greater ranges than current weapons systems. Several 
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laser weapons systems are currently in various stages of research and development. The 
airborne laser program (ABL), and tactical high energy laser (THEL) systems are 
examples of current projects designed for airborne and ground based applications. 
Applying this technology would greatly enhance the Navy’s abUity to defend against high 
end ASCM’s. A closer review of current ship missile defense systems will show that a 
directed energy weapon system will be a positive asset. 

1. Phalanx Close -In Weapons System (CIWS) 

The Phalanx close-in weapons system, shown in Figure 2 is a rapid-fire 20mm M- 
61A1 Gatling gun, which is capable of firing up to 4500 rounds/min of depleted uranium 
or tungsten armor piercing penetrators at a velocity of 1030 m/s [3a.]. The CfWS is a 
self-contained unit that can search, detect, track, and engage targets, but is typically 
integrated into a ships combat systems suite for additional fire control capabihty. 



Figure 2. Phalanx Weapons System (from [10]) 

The Phalanx weapon system is effective only at extremely short-range 
engagements, typically on the order of 1000 meters or less. The weapon can be fired for 
a li mited amount of time before the barrels begin to overheat; so that the engagement 
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time for an incoming missile is lim ited to about five seconds or less. At the firing rate of 
4500 shots/min, only 375 buUets are fired at the oncoming missile. Additionally, the 
magazine only holds 1550 rounds and requires several hours of time out of service to 
reload [3a.]. Another problem associated with CIWS is its availabUity. Average 
operational availabihty of CIWS averaged 75% on Navy ships for FY 97, 98, and 99 [5]. 
That means that 25% of the time a ship would be without its last line of defense against 
an ASCM. While a small magazine and lim ited firing time are significant, the biggest 
disadvantage of the CIWS is that most of the shots fired will not hit the incoming missile 
due to the shght dispersion in bullet trajectories from the gun. This small dispersion 
angle translates to a large error down range; as a result, most of the shots miss the 
incoming target. Figure 3 shows a typical bullet dispersion pattern at a range of 1000 
meters. 
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Figure 3. Phalanx Bullet Dispersion at 1000 Meters 
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For the simulation shown, a total of 300 bullets were fred with a total of 5 hitting 
on a typical missile cross section of 20 centimeters. This equates to a 1.66 % hit 
probabUity at a range of only 1000 meters. The probabihty of hitting an incoming missile 
decreases rapidly as range increases. A plot of probabihty versus range is shown in 
Figure 4. 



0 500 1000 Range (m) 1500 2000 2500 


Figure 4. Probabihty of Hit vs Range for a Single Buhet 

Lethahty testing shows that destruction of a cmise misshe requires multiple hits. 
The first step in determining the number of hits taken by a misshe during its approach is 
to estimate the number of buhets that the misshe whl encounter in a 100-meter range of 
its flight path. Assuming a misshe velocity of 400 m/s and a 75 shot/sec firing rate for 
the gun, the bullets encountered in each 100 meter section of path is calculated to be 26 
buhets. 

Using the predicted hit probabhities versus range and the number of buhets 
encountered in 100 meter increments, the typical khl range can be determined. Assuming 
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it takes an average rf 6 hits to “kill” a missile, then a typical range at which a missile wiU 
be destroyed is about 400 meters. The destroyed missile may break up into several 
smaller fragments, each of which could have sufficient kinetic energy to strike the ship 
and cause significant damage to personnel and unprotected equipment such as radar and 
antennae. The CfWS is the current last hue of defense against ASCM’s employed by 
many US warships. Based on lim ited engagement time, small magazine capacity, low 
probability of killing a missile at ranges outside of 500 meters, and high probabihty of 
damage from missile debris, the CfWS abihty to defend a ship against technologically 
improved missiles is rapidly diminishing. A directed energy weapon would improve 
engagement range and remove the dispersion problems encountered by the Phalanx 
weapon system 

2. Rolling Airframe Missile (RAM) 

The RAM is a point defense weapon system designed to engage a threat missile at 
intercept ranges out to 9.6 km [3b.]. It utilizes a Mach 2 missile that combines the 
infrared seeking of the Stinger missile with the warhead, rocket motor, and fuse from the 
Sidewinder missile for ASCM defense [4]. The typical configuration found on ships is 
the MK 49 launcher system, which has 21 missiles in cells ready to launch. The system 
receives input from the ship’s fire control system. Once the target information is 
received, the launcher will line up at the appropriate direction and elevation and fire a 
missile. The missile uses an RF seeker to acquire the target, and when sufficient IR 
signal is received from the target, shifts to IR mode for the terminal phase of flight. The 
missile is able to rapidly maneuver (up to a max im um acceleration of 20 g’s) to counter 
evasive measures by the incoming target. 

The RAM has a high success rate against ASCM’s, but it is not 100%. Time of 
flight to intercept a target allows range to close rapidly. If a RAM fails to kill the target, 
there will be tittle time to get the next RAM launched and the target will close to a short 
range, endangering the ship. Therefore, a preferred tactic is to fire two RAM at an 
incoming target to increase kill probability, utilizing the shoot-shoot-look method rather 
than the shoot-look-shoot method. This strategy gives a better probability of missile kill, 
but causes a rapid depletion of available weapons munitions and is expensive. The RAM 

system offers a greater engagement range and a higher hit probability than the CIWS. 
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Both systems are often used together such that the RAM provides longer range coverage 
and CIWS is used for short range defense against any missiles that get past RAM. 

3. Airborne Laser System (ABL) 

The ABL is an example of a directed energy weapon being designed to defend our 
forces against theater baUistic missile attacks. The ABL will engage a ballistic missile 
during the initial boost phase of flight and destroy the missile directly over, or near the 
launch site. The system consists of a megawatt class chemical oxygen iodine laser 
(COIL), beam control, transport section, and beam director. The system is mounted on a 
modified 747-4()0F aircraft with the beam director located in the nose cone section. 

The COIL is a continuous beam laser that operates by injecting chlorine gas into a 
reservoir of hydrogen peroxide liquid to excite oxygen. Iodine gas is then injected into 
the excited oxygen to produce excited iodine that decays to its ground state by emitting a 
photon of wavelength 1.315 |im [17]. The photons are amphfied in the lasing cavity and 
the energy is dehvered through the beam transport system to the beam director. The beam 
director contains adaptive optics that use piezoelectric actuators that adjust the beam 
shape to minimize the affects of atmospheric distortion. The platform will operate around 
45,000 feet to help minimize the atmospheric absorption and turbulence, to give the ABL 
an effective range of over 400 km [18]. 

The ABL system is evolving rapidly, and significant advances in laser technology 
have reduced the weight of the laser from 5500 lbs to just over 3000 lbs [4], not including 
the weight of beam control, beam director, or tracking systems. The program currently 
has one modified 747-400F that has the turret nose cone installed, with laser installation 
in progress. An expected operational test of the ABL is scheduled for sometime late in 
2003, with an engagement of a Scud-like missile [17]. 

Advances in technology do not come cheaply. The estimated total program cost 
for the ABL, to include the first two aircraft, is over six billion dollars with each 
additional aircraft costing in excess of $500 mi lli on [18]. In addition to the system cost, 
the ABL has an estimated cost per engagement in the range of several thousand dollars a 
shot. The ABL does have lim itations on the amount of chemical fuel that can be carried 
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but current estimates are that the system will be capable of firing the weapon a minimum 
of 20 times, with that number likely to grow in the future. 

This program shows that directed energy weapons are becoming a reahty. There 
are several weapons systems in various stages of development. Some other examples are 
the tactical high energy laser (THEL), and mid-infrared advanced chemical laser 
(MIRACL). High energy laser (HEL) weapons are likely to become a key component in 
combat systems over the next few decades. These new weapons provide speed of light 
capability with reduced time to kiU and confirm that a target is destroyed. The abihty to 
engage and kiU anti-ship missiles at a greater range than current systems is a realistic goal 
with the use of a directed energy weapons. Directed energy weapons will provide greater 
protection for the ship, and eliminate collateral damage to the ship from missile debris 
striking the ship. One type of directed energy weapon being considered for missile 
defense is the free electron laser (EEE). 

B. BASIC DESCRIPTION OF A FREE ELECTRON LASER (EEL) 

1. History of the Free Electron Laser 

The concept for the free electron laser was first introduced by John Madey in 
1971 [6], and in 1976 he demonstrated the first successful experiment using an EEL 
amphfier [7]. Interest in development of the EEE increased rapidly and by the 1980’s 
research was heavily funded by the Strategic Defense Initiative. Development of EEE 
designs and technology have continued to the point where high power EEE ’s are now 
being considered for weapons apphcations again, including anti ship missile defense. 
Thomas Jefferson National Accelerator Eacihty (TJNAE) in Newport, Virginia has an 
operational EEE that has demonstrated an average power output of over IkW and is 
currently being upgraded to 10 kW. Euture plans include modifying the laser to achieve 
an average output power of 100 kW with a projected demonstration in 2005 [9]. The 
technology for the EEE has advanced rapidly to the point where a shipboard EEE missile 
defense system is possible within the next decade. 

2. Electron Laser Basic Description 

The free electron laser operates by passing a relativistic electron beam through an 
undulator that produces a periodic magnetic field developed by a series of opposing 
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magnets. The spatially oscihatmg magnetic field causes the electrons to foUow sinusoidal 
trajectories. The acceleration of the electrons in the transverse direction causes them to 
emit radiation along he axis of the undulator. The fight is contained between two mirrors 
that make up the optical resonator cavity. One mirror is partially transmissive to allow for 
fight to escape the optical cavity and be used for the designed application of the laser. 
Figure 5 shows a basic FEL cavity. 



Figure 5. Oscillator FEL Cavity 


As the fight is reflected between the mirrors, it sets up an optical mode that is 
traveling in the same direction as the electron beam. The electron beam is not a 
continuous beam, but rather a series of pulses separated by a specific distance. This 

separation allows an electron pulse to enter the undulator in synchronism with the 

reflecting fight pulse. As the electron pulse travels through the undulator in the presence 
of the fight pulse, the electric field of the fight interacts with the electrons, causing 
stimulated emission of more fight. Over many passes, this creates a coherent fight beam 
that continues to grow as the fight reflects between the mirrors and interacts with more 
electron bunches. As the light intensifies in the cavity, saturation is reached where the 

output of the light is equal to the fight created in the cavity. 

The undulator and optical resonator cavity are where the laser beam is produced, 
but they are only a small part of the system that makes up the complete FEL system. 
Figure 6 shows a typical FEL system in a ring configuration with energy recovery. 
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FREE ELECTRON LASER 



Figure 6. PEL Ring Configuration with Energy Recovery 


The electron beam path is shown in blue, except where it overlaps the laser beam 
path shown in red. The electrons are generated in the injector by photoemission or 
thermionic emission, depending on the type of injector used, and accelerated to 
approximately 10 MeV by a strong electric field set up at the photo injector gun. The 
electrons are produced in a series of short pulses, typically on the order of milli meters in 
length, at a pulse frequency that is typically in the lOO’s of MHz range. The electron 
pulse leaves the injector and enters the hnear accelerator where several hundred MeV of 
electron energy is gained. The accelerator uses radio frequency (RE) cavities to produce 
large electric fields that accelerate the electrons to their final energy, which can be as 
large as several hundred MeV. The electron beam is then directed into the undulator by a 
series of bending magnets. Inside the undulator, a few percent of the energy in the 
electron beam is converted into radiated energy in the form of laser light. The electron 
beam then leaves the undulator and is directed back into the linear accelerator. The 
electron beam enters the accelerator 180 degrees out of phase with the RE field in the 
accelerator and gives a large portion of its energy back to the cavity while the beam is 
decelerated by the RE field. The low energy electron beam exits the accelerator and is 
then directed to the beam dump, where the remaining energy of the electrons is 
dissipated. 
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This energy recovery technique has two positive affects on the system. First, it 
increases the overall efficiency by regaining energy from the electrons before dumping 
them. Second, it reduces the electron energy to less than 10 MeV, which is lower than 
the energy required for neutron generation; this prevents neutron irradiation from causing 
activation of materials in the beam dump and minimizes shielding requirements. 

3. Other FEL Configurations 

The proposed FEL for shipboard defense is an oscillator configuration with 
energy recovery. There are several other types of FEL configurations that may be used 
for other apphcations. The linear oscillator FEL with energy recovery, shown in Eigure 
7, removes the need for bending magnets. To accomphsh energy recovery, a second 
linear accelerator and beam injector are used to inject electrons in the opposite direction 
of the first accelerator. The hnear accelerators are positioned such that electrons coming 
out of the undulator are 180 degrees out of phase and give up energy into the rf cavities. 
With two opposing accelerators, the tight beam is amplified in both directions due to the 
interaction with electron bunches from both accelerators. The additional space required 
for a second linear accelerator make this configuration impractical for shipboard 
application. 



Injector Undulator Injector 


Optical Resonator Cavity -1 

Eigure 7. FEL Linear Configuration with Energy Recovery 

Eigure 8 shows an amplifier FEL configuration that is a single pass high gain 
design which extracts as much energy from the electron as possible during the pass. The 
disadvantage with the amplifier design is that due to electron beam dispersion from large 
energy removal, energy recovery is not possible. The higher energy electrons create 


13 





















higher radiation levels at the beam dump and will have sufficient energy to cause neutron 
production and subsequent irradiation of material in the beam dump. 


Injector 


Linear Accelerator 




-/ F^^Laser Be^ m 

Laser 

Cavity 


Figure 8. FEL Amplifier Configuration 


In addition to the configurations discussed above, there are many other possible 
designs. Similar systems without energy recovery are common, but advances in 
technology have made energy recovery more attractive based on the increased efficiency 
of the system. 

C. FEL FOR MISSILE DEFENSE 

The current missile defense systems aboard Navy ships are rapidly falling behind 
the developments in stealth and speed of anti-ship cmise missile technology. A high 
energy laser (HEL) system, such as a FEE, would be a positive asset for many reasons. 
The dehvery of energy at the speed of tight to a target reduces the range that an incoming 
threat missile can close on a ship before being destroyed. Eigure 1 showed that the 
farther away a missile is destroyed, the less probable it is that a missile fragment will 
strike the ship. In addition to reduction in time of flight to the target, a HEE system 
allows for instant verification of a target kill and then allows for rapid reassignment to a 
different target if necessary. There are many types of HEE systems that are currently 
being developed for applications in the battles spaces of the future, but the FEE offers 
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many advantages for shipboard applications that other HEL systems are not capable of 
matching. The EEL also offers advantages over conventional type weapons such as 
RAM, and CIWS. 

1. Advantages of the FEE 

a. Supply 

Euture plans for all-electric drive ships will make large amounts of 
electrical power available for weapons apphcation. With this electrical energy available 
for use by the weapons system, a high power EEE would be able to operate and engage 
targets for as long as the ship’s electrical system can supply power. Conventional anti¬ 
ship missile defense systems that rely on an inventory of missiles or ammunition, and 
chemical laser systems that require a chemical source, face the problem of miming out of 
the required munitions and require an out of service period for reloading. The EEE 
would be never have to be removed from service to undergo reloading operations, 
making it much more rehable. 

b. Tunable Wavelength 

The EEE operating wavelength is tunable, which gives the abihty to 
exploit wavelength-dependent atmospheric propagation windows. By minimizing 
atmospheric absorption, the effective range is increased and the time required to damage 
a target is decreased. Other types of lasers such as the chemical laser, eximer lasers, and 
sohd state lasers can only operate at only one specific wavelength. One example of a 
non-tunable laser is the MIRAGE system, a deuterium-fluoride laser that operates at 3.8 
pm [8]. The wavelength (X) is determined by the energy levels excited in the chemical 
reactions of the lasing process making it impossible to change the wavelength of the laser 
to minimize atmospheric absorption. 

c. Reliability 

The EEE design requires a minimum amount of maintenance since there 
are no moving parts. The components are all electrical, with the exception of some 
support equipment such as the refrigeration/coohng units and the beam director. The 
EEE at TINAE is currently undergoing an upgrade to 10 kW and is projected to have 
10,000 hours of continuous operation before a scheduled maintenance period. This type 
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of reliability wiU be a great improvement over current existing systems that require more 
frequent periodic maintenance and repairs. The Phalanx CIWS is, on average, only 
available 76 % of the time for a variety of reasons that include hydrauhc problems, lack 
of onboard repair and preventive maintenance parts, and overhauls [5]. In addition to 
system rehabihty, the PEL offers another distinct advantage over missile systems. A 
beam of hght has very httle chance of suffering a catastrophic failure before getting to its 
designated target, but a missile has many failure mechanisms that may prevent it from 
reaching the target. The PEL, with minimal maintenance requirements and high 
reliabUity, wiU greatly improve anti-ship missile defense capabihties. 

d. Exhaust 

The PEE converts energy from an electron beam into a high energy hght 
beam without the use of a chemical fuel. There are no noxious gases or plumes that must 
be vented off ship, u nlik e a chemical laser or missile system. This minimizes controls for 
hazardous materials and e lim inates the need for expensive clean up in areas that are 
exposed to hazardous chemicals. 

e. Mission Flexibility 

The PEE is an adaptable weapons system capable of taking on many roles. 
Unlike other defensive weapons systems like the CIWS and RAM that are specifically 
designed only for missile defense, the PEE would be capable of defending against many 
other threats such as small ships, aircraft, and jet ski’s. In addition to a defensive role, 
the PEE can also be used for precision strike capabihty in littoral waters. The ability to 
target and destroy key equipment and assets with no collateral damage make the PEE an 
attractive alternative to a missile strike. The precision strike capabihty would naturaUy be 
limited to hne of sight targets that are located close to the coasthne so toat they would be 
within the range of the PEE. But there are now studies that consider using relay mirrors 
to re-direct the laser beam to the battlefield 

/. Operating Cost 

The cost per engagement for a PEE is much less than that of other 
weapons systems. The only cost incurred for firing the weapon is the cost of the fuel 
required to generate the energy used by the PEE, less than $2.00. In comparison, the cost 
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of one RIM-116A Rolling Airframe Missile (RAM) is $273,000 and one RAM is 
$444,000 [4]. A detailed cost analysis for the PEL is presented in section C. 3. of this 
chapter. 

2. Disadvantages of the FEL 

The FEL has many advantages as a weapons system, but there are also some 
disadvantages compared to other weapons systems. The first disadvantage is the high 
initial cost for the system, cmdely estimated around $50 mi lli on. Cost analysis for the 
FEL is presented in section 3. 

The other disadvantages of the FEL are related to beam propagation through the 
atmosphere. As mentioned in the advantages section, the atmosphere absorbs energy as 
the hght beam travels through it. Although the FEL’s abihty to tune the wavelength can 
minimize atmospheric absorption, it can’t e li minate it. Atmospheric absorption can 
change drastically based on environmental conditions. Fog, raindrops, and suspended 
aerosols cause additional scattering and absorption. As the particle density increases, the 
amount of absorption and scattering increase, resulting in a much lower beam intensity at 
the target. Fog, rain, dust, smoke, and other suspended aerosols greatly reduce the 
effective range of the FEL and may reduce intensity to the point that the weapon may be 
ineffective. Fortunately, these same adverse conditions affect the sensors on the threat 
missile. 

Another beam related disadvantage is thermal blooming, which occurs when a 
section of the atmosphere absorbs a small amount of energy and begins to shghtly 
increase in temperature. The refractive index of air is a function of temperature, and as 
the temperature changes in the center of the column of air being heated by the laser beam, 
it begins to act as a diverging lens and causes the beam to diverge. This problem is made 
worse with higher atmospheric absorption and beam intensities. Thermal bl(X)ming 
mcxiels use a critical bl(X)ming time, Xc , that estimates the time it takes for a column of 
air to be heated to the point that thermal blooming occurs. By ensuring that the beam 
does not heat a section of air for longer than the critical blooming time, thermal blooming 
can be avoided. Cross winds and slewing rate of the beam help to prevent thermal 


17 



blooming by changing out the volume of air in the beam. This is known as ehannel 
clearing. Thermal blooming will be diseussed in more depth in a later ehapter. 

3. Cost Comparison of the FEL vs Other ASCM Defense Systems 

A eost analysis of the FEL against current ASCM defense systems shows that the 
FEL system is a feasible alternative that may offer long term eost savings. The eost 
analysis is broken into two different areas: initial weapon eost and eost per engagement. 

Initial eost of the RAM system is $17.2M and ineludes $7.9M for the launeher 
and $9.3M for 21 missiles at $444,000 each [3b.][4]. This cost only provides for one 
launeher on the ship and does not provide for any additional missiles other than those 
loaded into the launcher. The initial cost for the Phalanx CIWS is $3.8M and ineludes 
the eost for one mount and 5000 rounds of ammunition [3a.]. The current estimate of a 
MW elass FEL is $55 mi lli on [9]. While the initial eost is sigrdfieantly more than the 
in itial cost of the Phalanx or RAM systems, once the FEL is installed, there are no 
additional eosts for ammunition, missiles, or stoekpile and storage of these munitions. 

Cost per engagement is based on a single engagement of an ASCM. Based on 
weapons doctrine, the normal tactie is to launeh 2 RAM at an incoming threat missile. 
Total eost for two RAM is $888,000. For the CIWS, an average of three (3) seeonds 
firing time is assumed to kill the incoming threat missile, whieh equates to 225 rounds of 
ammunition. Total eost for 225 rounds is $13,500 [3a.]. The FEL requires 
approximately five seeonds lasing time to put required energy on target to destroy it. For 
a FEL effieieney of 10%, this requires 50 Ml of energy from the ships eleetrieal system. 
One gallon of fuel eontains 113.5 Ml of energy. Take as an example, an LM2500 turbine 
whieh eonverts fuel to meehanieal energy at a rate of 0.435 lbs fuel/HP-HR. For a 
generator with a 90 % eonversion effieieney from meehanieal to eleetrieal power, one 
pound of fuel will yield 5.54 MJ of energy. To get the required 50 Ml of energy for the 
FEL, 9 lb, or 1.2 gallons of fuel are eonsumed. Using an average price of $1.40 per 
gallon, the total eost per enga^ment for the FEL is only $1.68. This low eost per 
engagement gives the abihty to use the weapon for five fire testing and training at low 
costs. The other two systems are rarely used for training due to large operating eosts, so 
personnel training is aceomphshed by simulations or walk through of proeedures. Based 
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on the low operational cost for the FEL, savings over decades of use make the system 
cost effective and in hne with other weapons systems. Table 1 shows a summary of the 
cost comparison for the FEL, RAM and CIWS. 



FEL 

RAM 

CIWS 

Unit Cost 

$55M 

Launcher - $7.9M 

Mount - $3.2M 



Missiles - $9.3M 

Ammo (10,000) - $0.6M 



Total - $17.2M 

Total - $3.8M 

Cost per 

2 Seconds fuel 

2 Missiles 

225 rounds 

engagement 

$1.68 

$888,000 

$13,500 


Table 1. Cost Analysis Summary 


Although the FEL is more expensive than RAM and CIWS for initial installation, 
over the hfetime of the weapon system, the FEL’s low cost per engagement coupled with 
the advantages hsted above make the FEL a valuable asset that will greatly improve 
protection for ships against ASCM’s. 


A system hfetime total cost is presented in table 2 using an assumed 20 year hfe 
with an average of 10 launches per year. 



FEL 

RAM 

CIWS 

Initial Cost 

$55M 

$17.2M 

$3.8M 

+ 

-1- 

-1- 

-1- 

Ammunition 

200 * $2.00= 

200*$888,000= 

200*$13500= 

Cost 

$400 

$178M 

$2.7M 

Total Cost 

= $55M 

=$190M 

=$6.5M 


Table 2. Lifetime Cost Analysis 
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m. DETAILED DESCRIPTION OE A EEL 


The PEL at TJNAF is currently being upgraded to an average output of 10 kW. 
Additional modifications are expected to increase that power output to 100 kW by 2005. 
These rapid advances in design and technology make a PEL with an average power of 1 
MW appear possible in the future. Simulations and modehng for the 100 kW EEL design 
at TJNAP are presented in this paper. Pigure 9 is a diagram of the proposed 100 kW PEL 
system at TJNAP [14]. 



Pigure 9. 100 kW PEL System Diagram (after [14]) 


There are two major subsystems to consider for the PEL: the electron beam 
control subsystem and the optical beam control subsystem. There are also additional 
auxihary systems associated with the PEL such as refrigeration, coohng water, and 
shielding. 

A. ELECTRON BEAM CONTROL SUBSYSTEM 

The electron beam control subsystem is composed of the injector, linear 
accelerators, beam dump, and beam transport subsystem. The electron life cycle can be 
traced in Pigure 9. The electrons are generated in the injector, accelerated to high 
energies in the accelerator modules, pass through the undulator where they give up a 
small amount of their energy to create the coherent laser beam, pass back through the 
linear accelerators where they give up the majority of their energy, and are then 
dissipated in the beam dump. 
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1. Injector 

The proposed injector design for the 100 kW PEL at TJNAF is a superconducting, 
500 kV DC, photocathode electron gun with a GaAs cathode that is driven by a green 
laser [15]. Figure 10 is a cutaway view of the proposed injector. 
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Figure 10. Injector Cutaway (after [15]) 

The injector generates the electrons in short periodic pulses by photoelectrically 
exciting the electrons off the cathode using a pulsed green laser. Once the electrons are 
removed from the cathode, they are rapidly accelerated by a strong electric field set up by 
the high voltage (HV) column. It is important to rapidly accelerate the beam as it exits 
the cathode to min im ize space charge effects due to Coulomb forces. The space charge 
forces cause the electron beam quality to degrade by causing the beam to spread in the 
transverse and longitudinal directions. The beam leaves the HV column and enters a 
series of three accelerator cavities that raise the energy of the electrons to greater than 7 
MeV before exiting the injector. 


22 




The gun is designed to operate with the following parameters to provide a high 
quality electron beam input to the superconducting RF (SRF) accelerators [12]: 


Pulse Repetition Rate 

a 

750 Mhz 

Peak Current 

Ipk 

270 A 

DC Voltage 

Vgun 

500 kV 

Electron Energy 

F 

J^inj 

>7 MeV 

Energy Spread 

AE 

< 1.0 % rms 

Transverse Emittance (normahzed) 


< 3 7C mm-mrad 

Electron pulse Length 

le 

0.1mm 

2. Accelerator 




After leaving the injector, the electron beam enters the hnear accelerators and is 
accelerated to relativistic energies. At such high energy, the space charge forces become 
insignificant and can be neglected in later calculations. The electron beam enters the 
accelerator at energies of approximately 7 MeV and is accelerated to a final beam energy 
of 210 MeV. The electrons enter the accelerator in phase with the RF field and are 
accelerated as they pass through the cavities. After leaving the accelerator the beam is 
directed to the undulator, where it gives up a small amount of its energy to the optical 
beam. After passing through the undulator, the electrons are returned to the accelerator 
180 degrees out of phase with the RF field, causing the electrons to decelerate by giving 
energy back to the RF field. 

The design parameters for the SRF accelerators are: 


Pulse Repetition Rate 

a 

750 MHz 

Acceleration Gradient 

(AE/L)acc 

20 MeV/m 

Number of Modules 


3 

Output Energy 

Eb 

210MeV 

Accelerator Efficiency 


60% 

Beam Dump 




The decelerated electron beam leaves the accelerators with less than 10 MeV of 
energy and enters the beam dump where the electrons and their energy are dissipated. By 
recovering the electron beam energy, the overall efficiency of the FEL is greatly 
increased. Due to degradation of beam quality, the electron beam can not be used for 
more than one pass through the undulator. 
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The beam dump consists of a block of metal into which the electron beam is 
dissipated, and a coohng system that removes the generated heat from the metal to a heat 
sink. Radiation shielding of the beam dump is required due to the production of gamma 
rays from Bremsstrahlung (braking radiation) produced as the energy of the electron 
beam is dissipated [13]. One advantage of the energy recovery system is that the final 
energy of the electron beam has been reduced below the energy required for neutron 
generation. The low energy gammas require less shielding than would be required for 
neutrons, which would activate adjacent materials and create additional radiation. 

4. Beam Transport Subsystem 

The beam transport system consists of piping and bending nagnets. The electron 
beam must be contained in high vacuum (10'^* Torr) to prevent interaction with gas 
particles. The beam is guided by a series of bending magnets that force the electron 
beam to change directions due to field interactions. The acceleration of the electrons due 
to these directional changes causes synchrotron radiation projected into a narrow forward 
radiation cone much like a fiashhght [13]. 

B. LIGHT BEAM CONTROL SUBSYSTEM 

Laser hght from the PEL is transported to a beam director cn the deck of the ship 
that sends the beam to the intended target. The major components of the hght beam 
control subsystem are the undulator, optical resonator cavity , and transport subsystem. 

1. Undulator 

The undulator, also caUed the wiggler, is contained in the optical cavity and is the 
heart of the PEL. Eigure 11 shows a detailed drawing of the optical resonator cavity and 
undulator [16]. 
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Figure 11. Optical Resonator Cavity and Undulator of FEL (after [16]) 

The electron beam is directed into the undulator where an oscillating magnetic 
field is set up using a series of alternating permanent magnets. The parameters of the 
undulator are: 


Undulator Length 

L 

2.88 m 

Undulator Wavelength 

Xo 

0.08 m 

Number of Undulator Periods 

N 

36 

Undulator Parameter 

K 

1.7 


2. Optical Resonator Cavity 

The optical resonator cavity is comprised of two mirrors that are spaced a distance 
apart so that the reflecting optical pulses interact in phase with the sequence of pulses 
from the electron beam. As the electron pulses travel through the undulator in the 
presence of a tight pulse, stimulated emission of tight occurs. The resulting coherent 
laser beam continues to be amplified by each successive electron pulse passing through 
the undulator. 

The mirrors are a key component in the FEL design and must be capable of 
withstanding high power densities. Based on an optical coupling output of 20 %, the 
power incident on the mirrors will be 500 kW to generate the 100 kW output. Power 
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density at the mirrors is well within lim its for this power, but could become a factor when 
increasing power output of the laser. Use of short Rayleigh lengths to increase mirror 
spot size, and therefore minimize power densities at the mirrors, is discussed in section 
V. One of the mirrors is partially transmissive to allow a percentage the hght beam to 
escape and become the output of the laser. The length of the optical resonator cavity (5) 
is 32 m. 

3. Optical Transport Subsystem 

The optical transport subsystem is similar to the beam transport subsystem. It 
consists of pipes and mirrors that transport the hght beam from the PEL to the beam 
director for use. In a weapon system, the hght beam will be guided to a beam director 
that wih aim and send the hght beam to the intended target. 

C. AUXILIARY SYSTEMS 

The FEE requires several support systems and an adequate power supply for 
continuous operation. Some of the required systems are the refrigeration system, fresh 
water coohng, shielding, and power supply. A brief description and purpose of each 
system fohows. 

1. Refrigeration System 

The refrigeration system is required to provide hquid hehum to the injector and 
accelerator. Coohng of the accelerator cavities to 2 K is required to eliminate resistance 
losses that would occur in the cavity wahs with such high electric fields. Supercooling 
the cavities makes them superconducting, thereby eliminating I R losses. 

2. Fresh Water Cooling System 

The fresh water coohng system is required to remove excess heat generated in the 
beam dump. The energy deposited in the beam dump must be removed to prevent 
overheating and damage to the beam dump. 

3. Shielding 

Operation of the FEE presents a radiation hazard that requires shedding. 
Bremsstrahlung radiation is generated in the beam dump and synchrotron radiation is 
generated in the areas where the electron beam direction is changed. The gamma rays 
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generated during FEL operation must be shielded to min im ize radiation levels in the 
general area. 

4. Vibration Control 

Vibrations due to couphng with external sources can cause mirror vibrations and 
motion of precisely ahgned equipment. To reduce the effects on the system, vibration 
isolation mounting is used. The FEL resonator mirrors require very precise ahgnment 
and must be maintained within a few microns for proper FEL operation. An active mirror 
ahgnment system is required to maintain ahgnment of the laser cavity and ophcal 
transport mirrors. 
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IV. FEL THEORY 


The free electron laser uses a relativistic electron beam, an undulator, and an 
optical resonator cavity to produce a high power, coherent, laser beam. Early theoretical 
descriptions of the FEL utilized quantum electrodynamics, but the development of a 
classical approach later proved to be both accurate and easy to understand [19]. The 
classical FEL theory will be described in the following sections. 

A. RESONANCE 

In order for an FEL to have gain, a net transfer of energy from the electron beam 
to the optical wave must occur. The optimum energy exchange between an optical field 
and electron occurs when one wavelength of light passes the electron in one undulator 
period. This is called the “resonance condition” and can be demonstrated as a race down 
the undulator between a photon and an electron where the photon wins the race by one 
optical wavelength. Figure 12 shows a diagram of the electron-photon race [21]. The 
optical wavelength X is shown in blue, the electron is red, and the undulator wavelength 
?io is in green. 



Figure 12. Electron - Photon Race 


The relationship between the undulator wavelength and the optical wavelength 
can be developed from the resonance condition. For a given reference electron, a 
wavelength of fight will travel a distance X ahead of the electron over one undulator 
period. Take At to be the time it takes an electron to travel through one undulator 
wavelength Xo at speed v^. This is equal to the time required for the photon to go a 
distance Xo-i-X at the speed of light c, so that 
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( 4 . 1 ) 


V, c 


Solving for X with = v/c gives the resonance condition, 


I 


^1-p ^ 

vX y 


Ko- 


The Lorentz factor is 


(4.2) 


1 _ 1 


(4.3) 


where K= e Ei, Xo/(2 n me c^) is the dimensionless undulator parameter and = KA{. 
(to be derived later in the electron dynamics section). Inserting K/^ for P_l into equation 
4.3 and solving for P^ we get 





1 - 


l+K^ 

2y" 


(4.4) 


The last approximation is made using the binomial expansion since K is of order 
unity and » 1. Inserting equation 4.4 into equation 4.2, the optical wavelength 
becomes 


1-1 + 


X=- 


l + K^ 

2j^ 


V 


7 


1- 


V 

[ 2y^ JJ 


K- 


(4.5) 


Since l + K^/ly^ « 1, the denominator is ~ 1 and the resonance condition can be 
written as 


1 = 




ir 


Xn . 


(4.6) 
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Equation 4.6 demonstrates that for a given undulator, the wavelength of hght 
becomes a function of the electron energy Eg = y me and the undulator magnetic field 
through K. 


B. ELECTRON MOTION 

1. Spontaneous Emission 

The electron trajectories in the undulator determine the character of EEL 
interactions [19]. Initially, there is no optical field in the resonator cavity; it is created by 
spontaneous emission from the electrons traveling through the undulator. Spontaneous 
emission of the electron occurs when the electron is accelerated - in this case the 
electrons oscillate transversely due to interaction with the magnetic field of the undulator. 
The power radiated by the relativistic electron, Pg , is given by the relativistic Larmor 
formula [21]: 


3 c 


(4.7) 


where |3^= Kkoc/y is the transverse acceleration, ko = 27tAo, and e is the charge of an 

electron [20]. The photons emitted by the electron will have an energy E = hc/k where h 
is Planck’s constant. By using the resonance condition for k, ko = 2n/ko , and equation 
4.7, the energy of the emitted photon can be written as 


E 


ph 


Itick^y^ 

1 + 


(4.8) 


The number of photons emitted spontaneously by one electron in one pass 
through the undulator. We , can be determined from the power emitted by the electron, the 
energy per photon emitted, and the total time the photon is in the undulator 


W = 


P At 

e 

E 


(4.9) 


W = 


liiN aK\\+K^) 


(4.10) 
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where a = jtie ~ 1/137, is the fine structure constant. For a typical FEL with FhlO^ 
and K~l, an electron emits approximately one photon per pass in the undulator into a 
narrow forward cone of sohd angle However, not all of the photons are emitted into 
the coherent optical mode. Only a small percentage of the photons emitted into the 
cone are within the smaller solid angle, l/N'^ determined by the coherent optical mode 
[19]. The number of photons emitted into the coherent optical mode during one pass 
through the undulator is then 

W^^aK^(\+K^) . (4.11) 

The total photon energy emitted by the electron is many orders of magnitude less 
than the total energy of the electron, therefore, the momentum recoil has a negligible 
effect on the electron’s path. The process of spontaneous emission from the electrons will 
not generate large power outputs, but it is necessary for the start-up mechanism in the 
FEL. In a typical EEL, the number of electrons in a pulse is on the order of 10^°. The 
number of photons emitted into the coherent optical mode during the first pulse through 
the undulator is sufficient to estabhsh a classical optical field that is amphfied through 
stimulated emission of subsequent electron pulses. 

2. Electron Dynamics - the Pendulum Equation 

Now that an optical field has been developed in the FEL, the interactions of the 
electrons with the alternating static magnetic field of the undulator, , and the moving 

electric and magnetic fields of the optical wave, Bs and Es, can be analyzed. The fields 
for a hehcal undulator can be written as [21]: 


Bu = B^ (cos(/:(, z), sm(kf^ z), O) 

(4.12) 

Bs = (sin('P), cos('P), O) 

(4.13) 

Es= E^ (cos('P),-sin('P), 0) 

(4.14) 

in cgs units where 'P = (kz-(Ot+^) is the phase of the optical wave, 

wave number, and (|) is the initial optical phase at t = 0 and z = 0. 

k = (O/c is the optical 
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The forces acting on the electrons are given by the Lorentz force equation, and the 

equation for the change in electron energy Eg , 

F =i^ = -e[E, + ^xB] 

dt ^ ’ 

(4.15) 

—^ = -e V- Es . 
dt 

(4.16) 

Using the relations, p = y m^v=(3 y m^c , and Eg = y mgC^, 

be rewritten. 

equations 4.15 and 4.16 can 

^(e;+Pxb) 

(4.17) 


(4.18) 


dt c 

Substituting the fields into equation 4.17, the transverse components can be written as 


- —(Es.+ ^XBuy+^XBsy] 

dt c ^ ' 

^^^ = -^(£.(l-|3,)cosm-5 |3^sin(^„ z)) (4.19) 

dt c 

and 

dt m ' 

e 

^^^ = -^(-£,(l-13jsin('P)+5„(3,cos(^o z)). (4.20) 

dt me 

e 

For relativistic electrons, 1- (3z «1, equations 4.19 and 4.20 can be simplified, and 
Jfv B 1 -e 5 B 

written together as- — = - - —^(-sin(fco z), cos(fco z), 0) . (4.21) 

dt c 

Integrating equation 4.21 gives (assuming constants of integration are zero for perfect 
injection into the hehcal orbits) 
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(4.22) 


Pi =— ^ (cos(^o z), 5in(fco zXO). 

y c 2 % 

Defining the dimensionless undulator parameter as K= e Bu Xo /(2 n ntg c^), equation 4.22 
becomes 

P = —(cos(A:o z), ^in(A:o z), 0). (4.23) 

Y 

The transverse motion of electrons in the helical undulator is described by equation 4.23. 

Substituting the electric field from equation 4.14 into equation 4.18, the second Lorentz 
equation becomes 

^ = (P.,P,,Pj-£.(cos'T,-sin'T,0). (4.24) 

at m„ c 


Inserting equation 4.23 for into equation 4.24 and using the trigonometric identity 
cos(a+h)=cos(a) cos(h) - sin(a) sin(h), equation 4.24 becomes 


dt 


cos(^oZ + ^) 

y c 


dy . e K E . . 

— =y = -^ cos(^ +(|)) (4.25) 

dt y c 

where ^ = (k+ko) z - CO t is the electron phase. 

By rearranging equation 4.4 and taking the time derivative, a relationship between 
y and P^ can be developed; 

|(y-^(l-ii:^) = l-p:) 

-2y-^y (l+ii:^) = -2 p^ p^ 


Y^IIPA 

y (l+iY^)' 

Taking the first time derivative of the electron phase, we get 


(4.26) 
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^ ={k+k^)v^ -0) = (A:+fco)c(3^ -0), 
and the second derivative is 

C=(^+^o)c|3,. 

Solving for (3^ yields 

(k+\^)c 

Substituting equation 4.29 into equation 4.26, 

Y_ 

Y {\+K^^[k + kQ) c 


(4.27) 


(4.28) 


(4.29) 


(4.30) 


Using the resonance condition (equation 4.6), (k+ko)c ~ kc=(0 since ko « k for 
relativistic electrons, and solving for Y , equation 4.30 becomes 

Y=^ (4.31) 

2o)o 

where (i)o=ko c = InAo- 

Equations 4.25 and 4.31 can be combined and solved for ^ 

. eKE^ 

Y=;:^ =-^ cos(C+4') 

2 (Oo y c 

C=2(0„^5^ cos(C+(^). (4.32) 

y c 

Equation 4.32 is the electron equation of motion and describes electron phase dynamics 
in the form of the pendulum equation. This equation can be written using the 
dimensionless time parameter x = ct/L, where t is the time and the undulator length is L. 
In a single pass through the undulator, x goes from 0 to 1. Eor clarity, derivatives with 

respect to dimensionless time will be indicated as dX fdx = X and the second derivative 
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will be X , where X can be any parameter. The pendulum equation can be rewritten as a 
second-order differential equation with respect to x as 


r-9 0 

y c 


f l2\ 


cos(^ 


(4.33) 


y y y y y 

where the L /c term comes from the substitution of d t=(L/c) d x. The group of terms in 
front of the cosine can be written as one variable, defining the magnitude of the 
dimensionless optical field a = I a I After some simphfications, it can be shown that 
the field is 


AnNe KE^L 

2 2 

y c 

so that the pendulum equation can be rewritten as 

^ =v = |a|cos(^ -l-(|)). 

where the dimensionless electron phase velocity v is equal to 



(4.34) 


(4.35) 


v = C =L[{k+k,)^^-k]. (4.36) 

Equation 4.36 governs the phase-space motion of electrons travehng through the 
undulator interacting with an optical wave of magnitude \a\. The cosine term determines 

if the electron gains or loses energy to the optical field. The electrons with phases from 
-7l/2to7l/2 gain energy from the optical field while, those with phases from 

71 / 2 to 37t/2 transfer energy to the optical field. As the electrons in the pulse interact 
with the optical field, they bunch over the optical wavelength. The optical field strength 
determines the rate of electron bunching. For |a|< 7l, the field is weak and the bunching 

is small, while for |a| > 7l, the field is strong and bunching is strong. 
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C. OPTICAL WAVE EQUATION 

The pendulum equation describes how electron motion is affected by the presence 
of an optical field. This section will develop the equations to show how the optical field 
is affected by the electron beam. The wave equation that governs the propagation of the 
optical wave is 


( 


V" 


V 


1 a"" 
7 ^ 


A = - 



(4.37) 


where A is the vector potential for a circularly polarized plane wave and is the 

transverse current density. For a slowly varying amphtude and phase of the optical beam, 
the vector potential can be written as a function of z and t. 


A{z,t) = -E^{ sin('P),cos('P),0) (4.38) 

(0 

and 


S =Vx A. 


(4.39) 


The transverse directions are neglected using the assumption that the electron beam 
radius is small compared to the optical mode. The spatial and time derivatives of 4.37 are 


d^A 2 dE.. 


dz^ k dz 


kAA 

dz 


|(cos('P),-sin(T0,0) +-E (-sin('P),-cos('P),0) 

k 


( 


kA± 

dz 


(4.40) 


and 


d^A 2 dE. 


dt^ k dt 


d(|) 

'dt 


- 0 ) 


(cos('P),-sin(TO,0) + -E (-sin('P),-cos('P),0) 

k 


d(|) 

'dt 


-CD 


(4.41) 


Using the assumption of a slowly varying optical phase and amphtude in time and space, 


dE, _ dE^ 

■« kE, , « k^ 


dz 


dt 


77 J A 

—«(oE , and —^ «C0(|) 
dz ' dt 
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the second order derivatives become negligible and are dropped from equations 4.40 and 
4.41. Using equations 4.40 and 4.41, equation 4.37 can be written as 


( 


2 A 




V 


A = 1 


dE, 1 dE 


dz 


- + - 


dt 


(cos('P),-sin('P),0) 


dz c dt 


(-cos('P),-sin('P),0) = 


471 


■7 . 


(4.42) 


Equation 4.42 can be simplified further by using a coordinate transformation z„= z + ct, 
and dimensionless time x = ct/L so that 


A+ii. 

dz c dt 


transforms to 


iA = A_ 

c dt Ldt 


making equation 4.42 become 


L dl 


(cos('E),-sin('I0,0)+2E 


L dl 


A 


471 — 

(-cos('E),-sin('IO,0) = -— J^. (4.43) 

c 


The transverse current density is the sum of the individual electron currents 


7i =-ecJ^P^5^(v-7j) 


(4.44) 


where 5^ (...) is the three dimensional Dirac delta-function and r. is the position of the f’’ 
electron. Inserting equation 4.24 for (3^ into equation 4.44, we get 

(cos(^ z), sm(ko z), 0)5^ (x - r.). (4.45) 

i y 

Equation 4.43 can now be decoupled into two separate equations, 
dE IzipIK 

-^ = - - —5^(v-/j)[(cos('E),-sin(TO,0) •(cos(fcoZ),sin(^oZ),0)] (4.46) 

and 
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IneKL y g 3 ^-cos('P),-sin('P), 0 )•(cos(fcoZ),sin(fcoZ), 0 ) 1 . ( 4 . 47 ) 

Upon completion of the dot product of the terms inside the brackets, the equations 
become 


dE, IneLK 


dx 


Y 


5^(^-r)cos(C +4)), 


(4.48) 


and 


d(|) _ -271 eKL 
dx yE^ 


Y 5^(^-q)sin(^+4)) 


(4.49) 


where ^ = (k+ko)z-(M. 


Assuming a constant density of electrons over the small volume element being 
considered, and using the average phase of the potential over the volume, the summation 
over i electrons can be written as 


Y S"(^-Y)cos(C+ 4 )) = p^(cos(C+ 4 ))), 


(4.50) 


and 


Y +4')= P.(sin(C+4')}. 


(4.51) 


where the average is denoted by (..) and Pe is the electron density. Using these relations, 
equations 4.48 and 4.49 become 

dE^ _ 2 n e L K 


dx 


Y 


-(cos(^ + 4 ))), 


(4.52) 


and 


^d(b In e K Lp, . ... 

Es^= -^(sm(^ +4))). 

dx Y 

These two equations can be formed from the real and imaginary parts of 


(4.53) 
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A(£_ e->) \, 

' Y ' ' 


(4.54) 


Multiplying both sides of equation 4.54 by 4neNKL/^mc^ gives 

A o _ 2 „2 XT 7^2 t2 


dx 


4 71 eN RLE 


1 2 

Y m^c 




8 71" e" Ni^" L" 


(e->. 


(4.55) 


By defining the dimensionless optical field, a, and the dimensionless current, y as 


a = a e 


II AneNKLE 8 7c" e" NiiT" L" 

\a\= -,-,- ^andj=— 


«,2 2 

Y m^c 


«,3 2 

Y 


equation 4.56 can be rewritten as the optical wave equation, 

a 


dx 


-a = a 




(4.56) 


The optical wave equation shows that the optical wave is dependent upon both current 
and average electron phase. With no current and/or no electron bunching, the optical 
wave will not change. 

D. GAIN 

The PEL oscillator exchanges energy from the electron beam to the optical beam 
over many passes, and can operate at low gain. Gain G is the fractional power change in 
the optical field per pass through the undulator. 


2 _ 2 
d] d'O 


(4.57) 


where ao is the optical field strength at the beginning of the undulator (x = 0)and ai is the 
optical field at the end of the undulator (c = 1). The electron beam must lose energy to 
the optical wave to achieve gain, thus one method to analyze gain is to determine the 
change in energy of the electrons. The energy of an electron is proportional to the 
electron phase velocity v given in equation 4.36 as 


v = ^ =(k + kQ)^f-(£) , 
or using dimensionless time x where dx =(c/L) dt, 


(4.36) 
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v = C =L[(k + k,)^^-k]. 


( 4 . 58 ) 


Using the approximation that k+ko ~ k (since ko«k), the change in the electron phase 
velocity is 

Av = L /t . (4.59) 


Using the resonance condition of equation 4.6, one can obtain the relation 


Av = 2 71 N 


( ^ 




By using equation 4.26 in a difference form, and solving for AP^ we get 


y (i+i^') 


AP.= 


AY(l+i^") 


(4.60) 


(4.61) 


Substituting into equation 4.60 and using the approximation that Pz ~ 1, we get a 
relationship between Av and Ay , 


Av = 4 71 



(4.62) 


The number of electrons in a small volume dV of an optical wave is given by 
Ne=PeFdV where F is the fillin g factor. The filhng factor is defined as the cross-sectional 
area of the electron beam divided by the cross-sectional area of the optical beam. Using 
equation 4.62, the average change of energy for an electron inside the undulator is 


Ay ~ 


y c"((v)-Vo) 
4 7C 


The energy contained in the volume dV is given by [21] 


de 


O —beam 


= ^dV. 
4k 


(4.63) 


(4.64) 
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Using the above relations, gain in the optical field is 


^ _ (9.FdV)l m, c’-({v)-v,)/(4 n N) 

lE^dvKin) -■ 

0-heam s I ^ ' 

Using the dimensionless parameter j, equation 4.65 can be simplified to 

G=^{vo-v)- (4.66) 

«o 

In weak fields, the changes in the electron phase, optical phase, and optical field 

are small so we can use the approximation that a ~ 0. Using the initial conditions 
v(0) = Vo, a(0) = ao, and applying perturbation theory, it can be shown that 


^ = ^0 +[sin(C 0 + VqX ) - sin( ^ 0 ) ] + 

Vn 


*^0 


1 


- - cos(2^o + 2voX) - cos(2^o) + cos(2voX) -1 - v^x sin(^ ^ )cos(Co + v^x) 


(4.67) 


+ ... 


To second order, the time average of the phase velocity is 

^2 

(v) = Vo +—^[2cos(VoX)-2+VoX sin(voX)]. (4.68) 

^0 

Substituting into equation 4.66, we arrive at the low gain equation 

G =-^(vo -v) = ^[-2cos(VoX ) + 2-VoX sin(voX)]. (4.69) 

ao Vo 

In the weak field, low current approximation, the gain is shown to be primarily a function 
of j and the initial electron phase velocity vo- Figure 13 shows a plot of the gain spectmm 
for the range of phase velocities from -12 to 12 in a weak optical field. The curve is anti¬ 
symmetric about Vo = 0 with a peak gain near 13% at an initial phase velocity vo ~ 2.6. 
The initial optical field is ao = 1, and dimensionless current is 7=1. 
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13 % 



Figure 13. Single Pass Gain (G) versus Initial Electron Phase Velocity (vo) 


It should be noted that at the resonance condition vo = 0, there is no net gain. The 
PEL must operate shghtly off resonance to have gain. 

E. PHASE SPACE 

The electron phase space is a plot of the electron’s phase velocity v versus the 
electrons phase Phase space evolution is helpful in showing the evolution of the 
electron as it moves through the undulator from x = 0 to x = 1. The ^ axis represents a 
section of the electron beam that is one optical wavelength long and is traveling at the 
resonance velocity vo- The plotted electrons move forward (or backward) based on their 
relative velocity in relation to the resonance velocity. 

Figure 14 shows a phase space plot for 20 sample electrons with an initial phase 
velocity vo = 0 and their evolution in phase space as they travel through the undulator. 
The electron position changes in color from yellow to red as he electron travels along the 
undulator. For this simulation, the electrons start out at the beginning of the undulator (c 
= 0) equally distributed in phase with an initial electron phase velocity vo = 0 (Yellow). 
The phase-space positions of the electrons at the end of the undulator (c = 1) are the 
final red dots for each electron. In this example, 10 sample electrons between -n/2 and 
nl2 increase their phase velocity (gain energy from the optical field) while the other 10 
electrons between 7t/2 and 3x72 decrease in their phase velocity (lose energy to the optical 
field). The net result is that there is no net gain G in optical field strength as shown in the 
upper right plot of Figure 14. It is clear that the electrons begin to bunch around ^ = nl2, 
which drives the optical phase, (|), shown in the lower right plot. 
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-ll/2 ^ 3x/2 0 T 1 

Figure 14. Phase Space Plot for a Low Gain FEL, vo = 0 


As shown in Figure 13, the peak gain for a low-gain FEL occurs when the initial 
electron phase velocity is vo ~ 2.6. A phase space plot with vo = 2.6 is shown in Eigure 
15. The phase space plot shows bunching of the electrons near ^ = 7t and shows that on 
average, the electrons have lowered in phase velocity thus giving a positive gain of 

G- 13%. 



^ Lj.i.Ui.i.i i 1 1 11 1 1111111 

-Z/2 ; 3Z/2 0 T 1 

Eigure 15. Phase Space Plot for a Low Gain PEL, vo = 2.6 
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F. 


fflGH CURRENT GAIN 


Equation 4.66 for a low gain EEL is only valid for weak fields |a|« 7C and low 
current j< n. In high currents j» n, the optical phase amphtude and phase change 
rapidly. Therefore the low current, low field approximation where \a\ is held constant 

and a = 0 can not be used for deriving the gain. The change in the optical field and gain 
in a high current EEL become exponential and are given by [19] 





a(x) 

1 ~ ^0 

1 2 J 2 

~ 3 



(4.70) 


and 


G(x) 


\ 

9 



(4.71) 


Equation 4.71 describes high gain at resonance as a function of dimensionless 
time X and current j, and is only valid for weak fields |a| << 7t with high current j »%. 
Eigure 16 shows typical high current EEL gain and optical phase spectra. 


**** Gain and Phase Curves **** 

j=200 a^=l N=20 

•' o 



Eigure 16. 





A phase space plot is useful in showing what happens in a high current FEL. 
Figure 17 shows a typical phase space evolution for a high gain FEL. An electron beam 
with an initial electron phase velocity Vq = 0 begins to bunch near the ^ = nil just as in a 
low current FEL. The bunching shifts the optical phase (|) resulting in a shift of the 
separatrix back in relationship to ^ as shown in the phase space plot. As a result, 
maximum optical amphtude growth occurs, and gain is exponential through the 
remainder of the undulator. Very high gains are possible in high current operations of the 
FEL. A gain of G~350 is shown in Figure 17 as a result of an initial current 7=200 in a 
weak field of ao=l. 


*** FEL Phase Space Evolution *** 

j=200 a =1 V =0 N=20 

•• o o 

Oq=0.5 i|=0.3% Av=10 



-121 _ 

-n/2 ^ 3n/2 Ox 1 


Figure 17. High Gain FEL Phase Space Evolution 


46 


















V. SIMULATION OF A SHORT RAYLEIGH LENGTH EEL 


A. INTRODUCTION 

In a low gain FEL, the cavity mirrors determine the fundamental optical mode, 
typically a Gaussian shape. At the narrowest part of the optical beam, the beam waist 
radius Wo is typically only a few millim eters and is usually located at the center of the 
undulator [23]. The beam radius will spread due to diffraction as it propagates in the z 
direction according to 



(5.1) 


where z is the displacement from the beam waist along the undulator axis, Zq = QW^/k is 
the characteristic spreading distance, called the Rayleigh length, and is the distance from 
the waist over which the beam waist area doubles. Transverse dimensions are normalized 
to ^JlXIti and longitudinal dimensions norm aliz ed to undulator length L. Thus, the 

normalized Rayleigh length is zo = Zq/L and the dimensionless waist w = WyjLk/% . 
Equation 5.1, written in dimensionless terms, where z = fx - 1/2)L, becomes 



The EEL is capable of producing extremely high average and peak power 
densities. The average electron beam intensities can be hundred’s of MW/cm and the 
intensity of the optical beam in the EEL can approach these intensities [22]. No mirror or 
optical material can withstand such high power densities without extensive damage. 
Therefore, the optical beam radius W must expand before reflecting on mirror surfaces to 
decrease the intensity on the cavity mirrors to acceptable levels. As shown by equations 
5.1 and 5.2, the spot size on the mirrors can be increased by either lengthening the optical 
cavity, and/or by shortening the Rayleigh length. Eor apphcations where space is limited, 
a short Rayleigh length is an attractive alternative to lengthening the cavity. 
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The TJNAF proposed 100 kW FEL wiU operate with an electron beam power of 
14 MW, thus requiring an extraction efficiency of 0.7% to reach the 100 kW output. 
Table 2 lists the proposed parameters for the 100 kW FEL. 


Parameter Symbol Value 


Electron Beam Eg 210 MeV 

Energy 


Pulse Repetition Q 750 MHz 

Rate 


Peak Current Ipk 270 A 


Electron Pulse k 0.1mm 

length 


Electron Beam rg 0.3 mm 

Radius 


Undulator Xo 8 cm 

wavelength 


Undulator Periods N 36 


Undulator length L 288 cm 


Undulator Parameter K 1.7 


Optical Wavelength X ~ 1pm 


Cavity Length S 32 m 


Resonator Quality Q 4.2 

Eactor 


Table 3. TJNAE 100 kW Parameters 

The resonator quality factor Q = 4.2 corresponds to a mirror transmission output 
of 1/Q ~20 %. Thus for a 100 kW output power there are 500 kW of optical power 
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impinging on each mirror. The power densities on the mirrors were calculated for the 
TJNAF PEL with dimensionless Rayleigh lengths of zq = 0.1, 0.2, 0.3, 0.4, and 0.5. 
Figure 18 shows the optical modes and mirror power densities on the mirrors for these 
different Rayleigh lengths. Reducing the Rayleigh length from ^ = 0.5 to ^ = 0.1 lowers 
the power density on the mirrors by 500 %. 



Figure 18. Mode Shapes and Mirror Intensities for Various Rayleigh Fengths 


In support of the proposed 100 kW FFF upgrade for TJNAF, numerous 
multimode simulations were mn to model and study the optical mode interaction with the 
electron beam. As the Rayleigh length changes, the optical mode shape changes. Due to 
resulting changes in the filling factor F, gain and steady-state power of the FFF may be 
affected. Using the pendulum equation to describe the electron motion and the optical 
wave equation to describe the optical field, three-dimensional {x, y, x) FFF simulations, 
were used to study these effects. Dimensionless parameters are used in the program to 
generalize results so they can be apphed to F FF ’s of various specific design parameters. 
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minimize numerical errors, and reduce equations by combining constants into meaningful 
variables. 

Figure 19 presents a three-dimensional simulation with a table of the 
dimensionless parameters for the simulation shown in the upper right hand block. 


FEL WAVEFRONTS 


8 


X 


-8 

8 

y 

-8 

20 


V 


-20 



Figure 19. Three-Dimensional Simulation Results For the TJNAF 100 kW FEL 


The dimensionless electron beam radius in the v: and y dimensions is = Oy = 0.4. 
The dimensionless betatron frequency is a)p= KkoLA{ =1 over the undulator length with 
the electron beam focused in the middle of the undulator at Xp = 0.5. Betatron motion 
describes the electron motion over many periods in the undulator with no tight present 
and is given as 

LB 

y(x) = yoCos((OpX )-i-^sin(a)px) (5.3) 

(Op 
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where yo and Gy are the initial position and injection angle of the electron as it enters the 
undulator relative to the undulator axis z.[23] The beam’s angular spread Oex = Oey = 0.16 
(rounded to 0.2 in figure) is determined using the matching requirement that = (0^ <3^ 
and is comparable with constant emittance results. The Rayleigh length for this 
simulation is zq = 0.2 

The plot in the upper left of Figure 19, |a(x,n)|, tracks the development of the 

optical mode over n=32 passes and shows how the optical mode develops in the cavity. 
The top center plot, |a(v,y)|, presents the wavefront cross-section as it exits the 

undulator at x = 1, and shows the electron beam (red) centered in the wavefront. The 
center plot, |a(v,x)|, is a cross-section of the optical mode in its final pass. The electron 

beam is shown in the undulator at each program iteration. In this simulation,to reduce 
computation time, the mirror separation was shortened to three times the undulator length 
instead of the actual separation of 11 times the undulator length. This does not change 
the result as the additional resonator length does not contribute to the optical field. The 
lower left plot, f(v,n), presents the electron phase velocity distribution and how it changes 
over 32 passes. The final electron phase-space plot is presented in the lower center plot 
and shows a final spread of Av = 24.5 which corresponds to an energy spread of Ay// = 
5.4% . The bottom right hand comer shows the development of gain G(n) and optical 
power P(n). Gain and power evolution are the parameters of interest for the simulations 
presented in the next two sections. 

B. WEAK FIELD GAIN SIMULATIONS 

Weak field gain simulations give insight into how a FEL wih start up from a cold 
cavity. Simulations were conducted for the proposed TJNAF FEL with variations of 
Rayleigh length from zo = 0. 1 to 0.5, electron beam radius from = Oy = 0.1 to 0.5, and 
initial electron phase velocity from vo = 1 to 15. Due to the large number of simulations 
required, the work was divided among several people, my contribution consisting of aU 
simulations for Rayleigh lengths zo = 0.1, 0.2 and analysis of data. Eigure 20 shows a 
plot of gain G versus initial electron phase velocity vo for the Rayleigh lengths of zo = 0.1, 
0.2, 0.3, 0.4, and 0.5 with an electron beam radius of Ox = Oy = 0.3. Peak gain for 
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Rayleigh lengths zo = 0.1 to 0.4 oceurred at an initial eleetron phase veloeity of vo = 4 and 
the peak for zo = 0.5 was at vo = 3. 



Figure 20. Weak Field Gain vs Initial Electron Phase Velocity 

The general trend of decreasing gain at higher phase velocities was broken for 
shorter Rayleigh lengths due to multiple optical modes that supported increased gain at 
higher phase velocities. Similar results were obtained for electron beam sizes Ox = Oy = 
0.1, 0.2, 0.4, and 0.5. 

Figure 21 shows a summary of weak field gain versus electron beam radius at the 
optimum electron phase velocity for Rayleigh lengths of zo = 0.1, 0.2, 0.3, 0.4, and 0.5. 
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Figure 21. 


Weak Field Gain vs Electron Beam Radius 


As the electron beam radius is increased, gain decreased in aU cases due to the 
increased number of electrons being outside the optical mode, reducing the gain. The 
combination of a small electron beam radius and short Rayleigh length increases the 
electron beam density causing amphfication of the optical mode and increasing the gain. 
Results from these simulations show that shorter Rayleigh lengths with small electron 
beam radius do not adversely affect weak field gain - in fact, they enhance it. 

C. STEADY STATE POWER SIMULATIONS 

The weak field gain simulations showed that the short Rayleigh length FEE has 
good gain and therefore will build optical power from the initial cold cavity with no 
initial optical field. The next step is to determine the final steady state power P 
achievable and the corresponding FEE efficiency p. Simulations for the proposed 
TJNAE 100 kW FEE were conducted by varying the Rayleigh length from zo = 0.1 to 0.5, 
the electron beam radius from Ox = Oj, = 0.1 to 0.5, and the initial electron phase velocity 
from Vo = 1 to 15 to determine final power and efficiency reached for each of the 
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combinations. These simulations use a strong initial optical amplitude to shorten 
program mn time. It has already been shown in the last section that weak field gains will 
build the optical field from cold cavity. Figure 22 is a plot of the FEL efficiency p versus 
initial electron phase velocity vo for five different Rayleigh lengths ranging from zo = 0.1 
to 0.5. 



0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 

initial Electron Phase Velocity 


Figure 22. FEL Efficiency vs Initial Electron Phase Velocity 


To understand how the electron beam radius affects efficiency, the highest peak 
power for each value of Rayleigh length (at the optimum initial phase velocity for that 
Rayleigh length) was plotted against the electron beam radius. Eigure 23 shows a plot of 
the FEL efficiency versus the electron beam radius. 
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Figure 23. FEL Efficiency vs Electron Beam Radius 

Eigure 23 shows that the FEL efficiency increases as the electron beam radius 
decreases. This increase is due to concentration of the electrons within the optical mode. 
Eigure 24 shows the maximum efficiency obtained for each Rayleigh length. 



Dimensionless Rayleigh Length zo 


Eigure 24. 


FEL Efficiency vs Rayleigh Length 
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In Figures 23 and 24, it is clear that the maximum efficiency for the FEL occurs at 
a Rayleigh length of zo = 0.3. This agrees with the predicted optimum fillin g factor 
occurring at Zq = l/(2'v/3) ~ 0.3 due to optimization of mode volume [19]. The 

maximum efficiency was 2% with a dimensionless electron beam radius of Ojc = = 0.1, 

and the initial electron phase velocity vo = 11 for zo = 0.3. 

With a short Rayleigh length and small electron beam radius, multimode 
oscillations were observed with steady-state power oscillating as much as 20%. Figure 
25 shows a short Rayleigh length simulation exhibiting multimode oscillations. 


|a(x,n) I 


-1 


40 


FEL WAVEFRONTS 


331 |a(x,y) | 23 



103 -8 


f (v,n) 


(C,V) 


j=360 

X =0 

o 

0)^=1.0 

Tg=0.5 

V =13 
o 


8 

|a(x,t) 1 

91 

y 



-8 

1 




ln(l+G(n)) 1.21 


-40 


n 


Figure 25. 


103 -K/2 ^ 3;r/2 0 

Multimode Oscillation Example 


103 


The oscillations are clearly visible as the modes compete resulting in an 
oscillating power at steady state conditions. The resulting efficiency, even with 
oscillations, is still greater than the required 0.7% for an output of 100 kW. 
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Based on the simulations, an PEL utilizing a short Rayleigh length provides 
sufficient gain and efficiency for 100 kW and greater operation, and allows for the power 
density at the mirrors to be reduced, thus preventing mirror damage. A short Rayleigh 
length design may be a step towards a compact high power PEL suitable for use as a 
shipboard weapon system. 


57 



THIS PAGE INTENTIONALLY LEET BLANK 


58 



VI. MIRROR VIBRATION AND FEL STABILITY SIMULATIONS 


A. INTRODUCTION 

In the previous section it was shown that short Rayleigh length FELs could 
operate with sufficient gain and efficiency to produce an output of 100 kW or greater for 
the proposed TJNAF FEL. The short Rayleigh length option appears to be a step in the 
right direction for increasing the power of a FEL to a MW class system. But operation of 
an FEL at short Rayleigh lengths introduces additional problems that need to be 
considered. The short Rayleigh lengths are possible only when the mirrors of the optical 
cavity are in a nearly concentric configuration. This arrangement may make the system 
very sensitive to mir ror vibration and misahgnment. A small angular change in the 
mirror position could possibly drive the optical mode to rotate such that it may reduce 
interaction with the electron beam reducing gain and power. 

To determine the effect of mirror vibration on the performance of the FEL, 
simulations were conducted for a proposed 1 MW FEL design. The parameters are 
similar to the 100 kW parameters presented in Chapter IV with the following changes. 
The dimensionless current is j = 210, corresponding to an electron beam energy of 185 
MeV with apeak current Ipk = 3.2 kA. The electron pulse remains 0.1 mm long with a 
pulse repetition rate of = 750 Mhz. The resulting power in the electron beam is 143 
MW, so that an extraction efficiency p =0.7% is required to achieve a 1 MW output. 
The dimensionless Rayleigh length is zo = 0.03 with an undulator length of L = 0.6 m and 
cavity length of 5 = 12 m. The optical wavelength for this design is X = 1 |im 

The simulations are conducted with the mirror at a fixed rotation Qsim- A static 
condition can be used for the mirror since the mirror vibration frequencies are on the 
order of kHz, but the fight only interacts with the mirror on the order of microseconds. 
Thus, the mirror appears to be stationary to the fight pulse. 

B. OPTICAL MODE TILT 

For a mirror that is tilted off the cavity’s optical axis by an angle 0, the optical 
mode will tilt by an amount tp. A corresponding shift in the mode spot on the mirror is 
given by 5y f5/2j(p. If the optical mode rotation is such that the mode no longer 
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contains the electron beam, the laser will not operate. Figure 26 shows a diagram of 
mode tilt and eorresponding spot shift at the mirrors. When the spot eenter shifts on the 
order of the mode radius (hy^w), the optieal mode will rotate tp =^max outside the eleetron 
beam and the laser will no longer funetion. 
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Figure 26. Optical Mode Tilt 


For and using the geometric estimate of 5}; = (Sll)ipmax and equation 5.1, it 

ean be shown that 
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With z=S/2 and Zo«S, equation 6.1 reduees to 


9 


max 



( 6 . 1 ) 


( 6 . 2 ) 


The optieal mode tilt as a funetion of mir ror tilt for a eold eavity (no eleetron 
beam present and thus no gain) is well known. Mirror tilt as a funetion of optieal mode 
tilt is given as [26] 

e„=(i+«)<p (6.3) 


where 0cc is the predieted mirror tilt for a cold cavity and the resonator parameter 
g=l-S/R where R = Sl2 + 2Z^/S is the mirror radius of eurvature. In the lim it that 
Zo«S, equation 6.3 beeomes 
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(6.4) 


Using the parameters for the proposed 1 MW design Zo = 0.018 m and S = 12 m, we see 
that dec = 18x10'^ tp. Thus, a small change in the mirror tilt aigle 0 should result in a 
large change in the optical mode tilt. 

C. SIMULATIONS 

The computer simulation used to model the mir ror tilt effect on PEL performance 
is the same three-dimensional PEL simulation used for the power and gain simulations 
presented in Chapter V. The same dimensionless parameters were used in these 
simulations and the program was modified to incorporate a tilt to one of the mirrors by a 
dimensionless amount 0m =d/(?i/nL/^^. 


The dimensionless mirror tilt angle 0m, and the initial electron phase velocity vo 
are varied in the simulations with the extraction efficiency p being recorded for each 
combination of 0m and vo. Eigure 27 shows a sample simulation output for a 
dimensionless mir ror tilt angle of 0m = 0.25, and initial electron phase velocity of vo = 10. 
The actual mirror tilt angle for the simulation is dsim = 200 prad with the corresponding 
optical mode tilt angle cpsim = 2300 prad. The output plots shown in Eigure 27 are the 
same types of plots that were described in Chapter V. The optical mode rotation can be 
seen clearly in the center plot. Note that the axes are not to scale since the y axis covers 
about 10 mm and the x axis covers 1.8 m. Thus, the optical mode rotation angle is 
greatly exaggerated. The optical mode rotation can also be seen in the top-middle plot 
where it is apparent that the electron beam is no longer centered in the optical mode 
wavefront at the end of the undulator x = 1. The extraction efficiency p for this extreme 
example of rotation was p = 1.3%, and is greater than the required 0.7% for a 1 MW 
output. 
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Figure 27. 


Mir ror Tilt Simulation Results 


D. RESULTS 

Figure 28 is a plot of efficiency p as a function of initial phase velocity vo for 
minor tilt angles 0^ = 0, 0.0825, 0.167, and 0.25 (0 = 0, 0.06, 0.12, and 0.18 mrad). 
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For small mirror tilt angles, efficiency goes up as the initial electron phase 
velocity vo increases until vg ~ 10 where there is a sharp drop in FEL efficiency. The 
drop in efficiency is due to a change in the optical mode shape going from a roughly 
Gaussian mode to a more complicated higher order mode stmcture. For larger mirror 
tilts, the drop is not as sharp and a more gradual change between modes was seen in the 
simulations. 
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Figure 29. Gaussian mode and Higher Order Optical Modes 



63 




















Figure 29 shows two cross sections for a mirror tilt of angle 0„, = 0.0825. The 
upper plot is for vo = 11.33 and shows a higher order optical mode at the mirrors with an 
extraction efficiency p = 0.47 % (less than the required 0.7% for a 1 MW output). In 
contrast, the lower plot for vo =1 shows a near Gaussian optical mode with a greater 
efficiency p = 1.06 %, than for the higher-order mode case. 

A comparison of the predicted maximum mirror tilt angle for a cold cavity Occ and 
the simulated maximum mirror tilt angle Bsim based on the same optical mode tilt cp 
showed that the cold cavity theory predicted a much smaller mirror tilt would produce a 
given optical mode tilt tpo (he 0cc«0i,m). For an optical mode tilt of cp = 2300 |irad, the 
cold cavity theory predicts a mirror tilt Occ = 0.04 prad , while our simulations actually 
gave Bsim = 180 prad. This shows that the optical mode is much less sensitive to mirror 
tilt than cold cavity theory predicts, indicating that the electron beam plays a dominant 
role in determining the final optical mode tilt. 



Figure 30 shows the FEL efficiency as a function of mir ror tilt. As expected, 
efficiency decreases as the mirror tilt increases. For very small tilt, almost no change in 


64 



efficiency was noted. As the tilt continued to increase, the efficiency began to decrease. 
Even for mirror tilts as large as 290 |irad (0m = 0.4), the efficiency was greater than the 
0.7% required for a 1 MW output. Current technology using active ahgnment techniques 
allows mirrors to be held stable to a tolerance of less than 0.1 |irad. Within this range 
there is clearly no noticeable change in EEL efficiency. Similar results were obtained in 
simulations conducted for the proposed 100 kW EEE at TJNAE [24]. Based on the 
results of the 1 MW and 100 kW simulations, mirror stabihty for a short Rayleigh length 
EEE does not appear to be an issue for operations when mirrors are stabi liz ed to a 
maximum tilt of less than 0.1 prad. 
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m POTENTIAL EEL EOR ASCM POINT DEEENSE 


A. INTRODUCTION 

The research conducted for this thesis has shown that a high-power FEL using a 
short Rayleigh length may be possible in the future. The need for an improved ASCM 
point defense weapon system was motivated in the first chapter and the advantages of the 
FEL over other conventional weapons systems and types of lasers were discussed. 

The FEL laser community is currently investigating high-power FEL’s for 
shipboard application and taking steps to raise achievable FEL power. The TJNAF is 
upgrading its existing FEL to 10 kW average output and will follow up with another 
upgrade to a system capable of 100 kW’s of average power. Several designs for scahng 
up to a MW class laser have been proposed and are currently being studied. There is stiU 
the question of how much power will be required b provide an adequate defense against 
ASCM’s? 

B. REQUIRED POWER OUTPUT FOR ASCM POINT DEFENSE 

Several factors must be taken into account to determine the required output power 
for a ASCM point defense weapon, these include (but are not li mited to) the energy 
required to destroy the target missile, acceptable dwell time on the target (length of time 
laser must remain on target to deliver sufficient energy to destroy it), and the losses due 
to beam propagation through the atmosphere. A logical way to determine the required 
laser power is to start at the target and work back to the ship. 

I. Energy Required to Destroy the Target Missile 

Determining the energy needed to destroy the target is not as easy as it sounds. 
There are several ways that the missile can be “destroyed”, by which it is meant that the 
missile or missile debris will not be able to strike the ship. Possible kill methods include 
warhead detonation, fuel tank detonation, stmctural damage, and sensor/guidance system 
interruption. Warhead detonation, fuel tank detonation, and stmctural damage are 
classified as “hard kills” because the missile is destroyed. Sensor and guidance system 
damage is classified as a “soft IdU” since the missile is stiU flying, but may be unable to 
dehver its payload. While a soft kill will generally require less energy, a hard kill will be 

used as the basis for determining the energy required to destroy the target missile. 
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There are many different kinds of threat missiles that the PEL system must be 
designed to destroy, and they all may have different areas of vulnerability. The warheads 
and fuel tanks are located in different areas for different types of missiles and there are 
different types of explosive, and fuel used. To simphfy the analysis, the required energy 
will be based on causing stmctural damage to the missile and allowing aerodynamic 
forces to destroy the missile. To cause sufficient stmctural damage so that the missile 
becomes unstable and breaks up due to aerodynamic forces, it is assumed that 0.5 hter of 
material must be melted. Many materials are used in the constmction of modem missiles, 
but aluminum is often used and will be assumed in these calculations. The energy 
required to melt the material is given by 

E ,=oV [ciT -T) + ^H ), (7.1) 

where Pm is the material density, is the volume melted, C is the specific heat, TJn is the 
melting temperature. To is the initial temperature, and AHm is the latent heat of melting. 
The values for aluminum are Pm = 2.7 g/crn^, C = 896 J/kg-K, Tm=855K, and AHm = 
4x10^ J/kg. Assuming an initial temperature of Tq = 400K, the energy required to 
destroy the missile is Emeit ~2. MJ. This calculation should be treated as an order of 
magnitude estimate. 

The estimated energy of 2 MJ agrees with experimental data collected in 

'y 

conjunction with the MIRACL program, showing that tens of kW/cm are required to 
destroy a missile with a dwell times of a few seconds [25]. For an average spot size on 
the missile of 100 cm , this gives energy in the MJ range. It will be assumed for the 

'y 

remainder of the calculations that about 2 MJ of energy deposited in a 100 cm spot size 
is sufficient to destroy a missile. This is by no means an exact number, but it gives a 
ballpark estimate for discussion. 

2. Propagation Losses 

As the laser beam travels through the atmosphere, it can lose energy in several 
different ways. The energy can be absorbed by the atmosphere, or get scattered in 
different directions, or, if the column of atmosphere that the beam is going through heats 
up enough, thermal blooming can occur. 
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a. Absorption and Scattering 

Absorption and scattering of the laser beam are related to atmospheric 
conditions and will be considered in this section. The amount of energy received at the 
target is a fraction of the energy that is transmitted from the ship and can be described by 
the equation Etarget = T Eshtp, where T is the atmospheric transmittance. The atmospheric 
transmittance T is given by 

T=e-^‘^ (7.2) 

where is the extinction coefficient and z is the distance the fight must propagate. The 
extinction coefficient is a sum of the absorption and scattering coefficients 

and is a function of wavelength. The scattering and absorption coefficients are dependent 
on the concentration and types of particles that are in the atmosphere and therefore vary 
greatly with different atmospheric conditions. The major contributors are CO 2 , ozone, 
water, and suspended aerosols. Figure 31 is a graph of the absorption and scattering 
coefficients of some of the major contributors to losses as a function of wavelength [26]. 
(Figure 31 uses the symbol o for the coefficients instead of a which is used in this 
paper.) 
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Figure 31. Extinction Coefficients vs Wavelength (from [27]) 


Figure 31 does not include the scattering and absorption due to water or 
CO 2 . When the effects due to water and CO 2 are included, it becomes apparent that some 
wavelengths propagate much better than others. Figure 32 includes the effects of CO 2 
and water and shows transmission windows where certain wavelengths have a much 
greater percent transmission than others [27]. 
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Figure 32. Percent Transmission vs Wavelength (from [27]) 


The abihty to take advantage of these atmospheric windows of transmission using the 
tunabUity of the FEL is a great advantage over other types of lasers. For any high power 
laser to successfully dehver lethal amounts of energy to a target, it must operate within 
one of these transmission peaks. Figures 31 and 32 are for a typical atmosphere over 
land. A maritime atmosphere contains a much larger concentration of aerosols than the 
atmosphere presented in figures 31 and 32. Figure 33 shows a plot of absorption, 
scattering, and overall extinction coefficient as a function of wavelength for a typical 
maritime atmosphere [28]. The lowest extinction coefficients occur at wavelengths of 
1.62, 1.25, and 1.06 pm. 


Maritime Atmosphere 


Marine Aerosol, Horizontal Path 



Figure 33. Coefficients of absorption, scattering, and extinction in a maritime atmosphere 

(from [28]) 
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The graphs presented in this section are for a typical atmosphere in clear 
weather at sea level. These graphs present only a small fraction of possible conditions 
and can not be used to satisfactorily describe the atmospheric absorption and scattering 
losses for all conditions in all areas of the world. Based on current available data for 
scattering and absorption, the best propagation windows in order of preference appears to 
be the 1.62, 1.25, and 1.06 |im wavelengths, but another factor must also be considered, 
thermal blooming. 

b. Thermal Blooming 

Thermal blooming occurs when the column of air through which the laser 
is propagating is heated causing the air to act as a diverging lens. The onset of thermal 
blooming is a function of the power density of the beam, absorption coefficient of the 
atmosphere, and the time that the laser beam is acting on the same column of air. For a 
given power density and atmosphere, the time until the onset of thermal blooming in a 
column of air is known as the critical blooming time Xc. The column of air through which 
the laser is propagating is continually changing due to cross winds and the slewing rate of 
the beam moving through the atmosphere to stay on target. One estimate used to 
calculate the critical blooming time shows that the critical blooming time is inversely 
proportional to the one third power of the product of the absorption extinction coefficient 
and beam intensity, Xc o«l/(aa I) • Figure 34 shows the peak transmission intensity as a 
function of wavelength including the effects of absorption, scattering, and thermal 
absorption for a crosswind of 10 m/s to account for channel clearing [28]. 

Maritime Atmosphere 



Wavelength (^m) 

Figure 34. Transmittable intensity through a maritime atmosphere (from [28]) 
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Figure 34 shows that the optimum transmission is at a wavelength of 1.06 
pm. Referring back to Figure 33, we find a value for the total extinction coefficient is ttg 
= 0.144 km * for 1.06 pm wavelength. Figure 33 assumes a “standard clear” visibility of 
23.5 km. Changes in the atmospheric visibihty will cause significant changes in the value 
of a. For visibihties of 60 km (exceptionally clear day) and 8 km (light haze) the 
extinction coefficients are 0.042 km* and 0.340 km * respectively. Table 3 shows the 
percentage of transmitted power from the ship (% Frans) that arrives at a target based on 
ranges of 1, 5, and 10 km for three different viewing conditions at the optimum 
wavelength of 1.06 pm. 


Viewing Conditions 




Exceptionally Clear 
(visibility = 60 km) 

95.9 

81.1 

65.7 

Standard Clear 
(visibihty = 23.5 km) 

86.6 

48.7 

23.7 

Light Haze 
(visibihty = 8 km) 

71.2 

18.3 

3.3 


Table 4. Percent Power Transmission as Function of Range at Various Viewing Conditions 

The table above demonstrates how atmospheric conditions can 
significantly affect the performance of the FEL. These predictions may radically change 
for different environments. Many maritime atmospheres contain different types and 
concentrations of aerosols making the propagation profile completely different than these 
models. The tunabihty of the FEL makes it ideal as a directed energy weapon in order to 
adjust to these propagation windows and dehver maximum energy to the target. 

C. REQUIRED LASER POWER OUTPUT 

To determine the required laser power output from the ship, several isseus must 
be considered. What kin d of dwell time is acceptable? What is the closest that I want a 
missile to get to the ship? How many missiles per engagement should be considered? 
The hst goes on. To determine the laser power at the ship, the following engagement 
scenario and assumptions were used: 
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1 . 


Engagement of two Mach 3.5 missiles with a detection range of 12 km and 
engagement beginning at a range of 10 km. 

2. Standard clear day, tte = 0.144 km'^ corresponding to a visibUity of 23.5 km. 

3. 2 MJ of energy “kills” the missile. 

4. 0.5 seconds required for retargeting after the first missile is destroyed. 

5. The last missile must be destroyed at 1000 m so that there is a less than 5% 
chance of fragments hitting the ship. 

For the scenario hsted above, the power of the laser at the ship would have to be a 
minimum of 1 MW. At 1 MW the first missile would be destroyed at 4.4 km with a total 
dwell time of 5 seconds to dehver the required 2 MJ of energy. The second missile will 
be destroyed at a range of 1km with a total dwell time 2.5 seconds. The difference in 
dwell times is due to the decreasing atmospheric losses as the missile gets closer. Figure 
35 shows a plot of the energy absorbed by each missile as a function of range for a 1 MW 
laser. 


Energy Absorbed vs Range 



Range(km) 


Figure 35. Energy Absorbed by Missile vs Range 

It should be noted that there are many assumptions used in making this final 
determination of the required power for a shipboard laser. This is only intended to show 
order of magnitude requirements for this apphcation, and indicates how a MW class laser 
is a probable solution for the growing ASCM threat. 
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D. 1 MW FEL 

The previous section showed that a 1 MW FEL could be a potential ASCM 
defense system. The system would be similar to the 100 kW FEL described in Chapter 
in including the electron beam control subsystem, hght beam control subsystem, and 
auxiliary systems. However, the size of the FEL will have to scale down while power 
will need to scale up from current FEL’s. There are several proposed MW class lasers 
being considered for the future. One such design was introduced in Chapter VI and 
simulations indicated a power output would be greater than 1 MW. 

For a weapon system design to be feasible on a Navy ship, it must meet certain 
volume and weight li mits. A comparison of the volumes and weights of three point 
defense weapons systems shows that the FEL is comparable in weight and volume to 
other systems. Table 4 shows a size and weight comparison for a FEL, Phalanx CIWS, 
and RAM system with a MK-49 launcher [2,3,4,9,10]. 


a 

FEL 

(2 Beam Directors) 

PHALANX 

(2 units) 

RAM system 

(2-MK49 launchers) 

Size 

FEL: 96 m^ 

Beam Director: 

16 m each 
Total: 128m^ 

57 m^ each 

Total: 114 m^ 

26 m each 

Total: 52 m^ 

Weight 

23,000 kg 

(estimate) 

6,170 kg each 

Total: 12, 340 kg 

Missiles: 42 @ 78.5 kg 

3300 kg 

Launcher: 2 @ 6117 kg 

12230 

Total: 15530 kg 


Table 5. Comparison of Size and Weight of Point Defense Systems 

Based on the advantages a FEL offers, the estimated size and weight comparisons, 
and the possibihty of being able to engage multiple ASCM’s successfully, the MW FEL 
is an excellent choice for an improved point defense weapon and shows promise as a 
weapon system for the future. 
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vni. CONCLUSIONS 


Rapid advances in free electron laser technology are making high power FEL’s 
possible. With the current upgrade of the TJNAF PEL to a 10 kW output scheduled for 
first hght in March 2003, and another upgrade to an output of 100 kW scheduled to begin 
in 2005, high power PEE’s are becoming a reality. To be a feasible weapons system 
capable of being deployed on Naval Combatant’s, size and weight restrictions must be 
adhered to. One of the limiting factors for sizing considerations is the power density on 
the optical mirrors. To maintain the power density on the mirrors at acceptable levels, 
the optical cavity must be long enough to allow for diffraction. Alternatively, a short 
Rayleigh length may be used to cause rapid spreading of the beam. Simulations 
presented in Chapter V showed that by using a short Rayleigh length PEE, power 
densities at the mirrors are significantly reduced with no significant reduction in the 
power output of the PEE. 

Por a short Rayleigh length PEE, the resonator cavity mirrors are in a nearly 
concentric configuration and are sensitive to misalignment and vibration. Small 
vibrations of the mirror may be sufficient to cause a decrease in the performance of the 
PEE, and if severe enough in magnitude may prevent operation of the PEE. Simulations 
presented in Chapter VI showed that as mirror tilt angle increased, the PEE efficiency 
decreased. However, it takes a mirror tilt several orders of magnitude greater than the 
active ahgnment tolerance of 0.1 |irad before the PEE efficiency is noticeably affected. 
Mir ror tilt, vibrations, and misahgnment within achievable tolerance limits wih not 
adversely affect the performance of a PEE. 

A MW class PEE appears possible in the near future. Designing a weapon system 
capable of providing quicker reaction time, and speed of hght delivery of lethal energy 
wih prove to be a valuable asset for the Navy, giving ships an improved ASCM defensive 
capabihty along with surgical strike capabihty in the httorals. Continuation of this 
research is needed to advance PEE technology to the point that a shipboard weapon 
system is possible. 
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